Nonperturbative renormalization group for the Kardar-Parisi-Zhang equation: 
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We present an analytical method, rooted in the nonperturbative renormalization group, that al- 
lows one to calculate the critical exponents and the correlation and response functions of the Kardar- 
Parisi-Zhang (KPZ) growth equation in all its different regimes, including the strong-coupling one. 
We analyze the symmetries of the KPZ problem and derive an approximation scheme that satisfies 
the linearly realized ones. We implement this scheme at the minimal order in the response field, and 
show that it yields a complete, qualitatively correct phase diagram in all dimensions, with reason- 
able values for the critical exponents in physical dimensions. We also compute in one dimension the 
full (momentum and frequency dependent) correlation function, and the associated universal scaling 
function. We find a very satisfactory quantitative agreement with the exact result from Prahofer 
and Spohn [l|. In particular, we obtain for the universal amplitude ratio go ~ 1.149(18), to be 
compared with the exact value go = 1.1504... (the Baik-Rain constant [2|). We emphasize that all 
these results, which can be systematically improved, are obtained with sole input the bare action 
and its symmetries, without further assumptions on the existence of scaling or on the form of the 
scaling function. 
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I. INTRODUCTION 

In their seminal work Q, Kardar, Parisi and Zhang 
proposed a stochastic continuum equation to describe 
surface growth through ballistic deposition which reads 

^^ = vV^h{t,x) + ^ {Vh{t,x)Y + #,f) (1) 

where rj is an uncorrelated white noise of strength D, 
{rj{t,x)r]{t',x')) = 2DS'^{x- x')Sit - t'), which mod- 
els randomness in deposition. The term v'V^h(t,x) pro- 
vides the interface with a smoothening mechanism. The 
insightful feature of Eq. ([T]) is to take into account the 
nonlinearity of the growth velocity through the inclusion 
of the term X{Wh{t,x))'^ which plays an essential role in 
the large scale properties of the height profile h{t, x). 

The KPZ equation ^ is maybe the simplest nonlin- 
ear Langevin equation showing non-trivial behavior [J], 
and as a consequence it arises in connection with an ex- 
tremely large class of nonequilibrium or disordered sys- 
tems [J| such as randomly stirred fiuid (Burgers equation) 
0, directed polymers in random media [6|], dissipative 
transport [3, Q or magnetic flux lines in superconduc- 
tors jj]. The KPZ equation has thus emerged as one of 
the fundamental theoretical models to investigate univer- 
sality classes in nonequilibrium scaling phenomena and 
phase transitions f4|. It is only recently, though, that a 
definitively-convincing experimental realization has been 
brought out for a one-dimensional interface, confirming 
detailed theoretical predictions |1CI| . 



The profile of the stationary interface is usually char- 
acterized by the two-point correlation function 



C(i, \x\) = {[hit, x) - /i(0, 0) - t{dth)f) 



(2) 



and, in particular, by its large scale properties. The 
KPZ growth leads to generic scaling. At long time 
T and large distance L, C assumes the scaling form 
C{t,L) oc t-^^/^ (/(L^^T~^^/^) without fine-tuning any 
parameter of the model. The scaling function g is uni- 
versal and has the asymptotics g{y) — > const, as y — >■ 
and g{y) ~ |y| as y — >■ cxd. x f-nd ^ ^^^ the univer- 
sal roughness and dynamical exponents respectively. In 
fact, these two exponents are not independent since the 
Galilean symmetry [Sj — the invariance of Eq. ([l} un- 
der an infinitesimal tilt of the interface — enforces the 
scaling relation 2: -I- x = 2 as long as A 7^ 0. 

The KPZ equation encompasses two distinct scenar- 
ios depending on the dimension of the interface. Above 
two dimensions, there exist two different regimes sepa- 
rated by a critical value Ac of the nonlinear coefficient 
d, Q. In the weak-coupling regime (A < Ac), the in- 
terface remains smooth, its properties are determined by 
the A = (Gaussian) fixed point — corresponding to 
the linear Edwards- Wilkinson equation [j| — with expo- 
nents X = (2 — d)/d and z — 2. In the strong-coupling 
regime (A > Ac), the nonlinearity becomes relevant and 
the interface roughens. In this regime, the exponents are 
not known exactly, and some important issues are still 
controversial, such as the existence of an upper critical 
dimension (for a recent discussion, see e.g. |ll|). This 



unsatisfactory situation has persisted up to very recently 
|12l | because the strong coupling phase of the KPZ equa- 
tion has remained out of reach of controlled analytical 
approaches. In particular, standard perturbation expan- 
sions [3, [ij, [ij] have been proved to fail at all order to 
find a strong coupling fixed point [15[. Some nonpertur- 
bative approaches have been devised, such as the mode- 
coupling (MC) approximation 7, 16 19], or the weak 
noise scheme [20|, but they are difhcult to improve in 
practice. (Some specific comments regarding the MC ap- 
proximation are put forward throughout this paper) . Let 
us finally mention the self-consistent expansion, which 
is performed on the Fokker-Planck equation and can be 
systematically improved (2ll - l23l |. It has yielded impor- 
tant results for the KPZ equation in arbitrary dimen- 
sions, where it predicts an upper critical dimension [11|. 
Note that it has also been applied to different extensions 
of the KPZ model, such as the non local KPZ equation 
2l| or the KPZ equation with colored noise [22]. 

In two dimensions and below, the situation is differ- 
ent, both from a physical and a theoretical point of 
view. Physically, the interface always roughens, there 
is no transition. On the theoretical side, exact results 
are available in one dimension [l|, y, [24l - [2a | . The criti- 
cal exponents x — 1/2 and z = 3/2 have been known 
for long since they are fixed by the existence of an inci- 
dental fluctuation-dissipation theorem in this dimension 
Q, but the scaling function g{y) (and other universal 
properties) have been computed exactly only very re- 
cently [l|, [j, |2^, |2g|. Note that there had been earlier 
attempts to determine the scaling function, in the frame- 
work of the mode-coupling approximation [I^ [23, [28J . 
In Ref. [231, a refined ansatz for the scaling function in 
one dimension is devised to solve self-consistently the MC 
equations and it turns out that the result compares quite 
accurately with the exact solution. The one-dimensional 
scaling function has also been obtained within the self- 
consistent expansion in one-dimension 23| and coincides 
in many respects with the MC result (see Section IVip . 
Regarding the MC theory, it is an 'adhoc' approximation 
which consists of resumming one-loop diagrams while dis- 
carding vertex corrections. The quality of the results is 
all the more surprising that it was shown that the con- 
tribution of the neglected terms are of the same order as 
those kept [lg| . The main drawbacks of the MC approach 
are that it strongly relies on the quality of the 'educated 
guess' for the ansatz - only available in one dimension 
up to now - and it cannot in practice be systematically 
improved by calculating higher orders. 

Recently, we have proposed an analytical approach to 
the KPZ equation based on nonperturbative renormal- 
ization group (NPRG) techniques [1^. This early work 
has shown that the NPRG flow equations embedded a 
strong-coupling fixed point in all dimensions, and it has 
thus provided a complete, qualitatively correct phase di- 
agram for the first time within a RG approach, as well 
as reasonable values for the strong-coupling fixed point 



In this paper we present a general and systematic 
framework for applying NPRG methods to the KPZ 
problem, strongly constrained by the symmetries of this 
model. We derive general Ward identities and introduce 
a 'covariantization' associated with the Galilean symme- 
try which, to the best of our knowledge, has never been 
reported before in the literature. Within this formal- 
ism, successive orders of approximation are easily made 
explicit. In the present contribution, we implement the 
minimal order in the response field of this approximation 
scheme. We review the (simplified) first account of it pre- 
sented in [l2| - postponing its complete revisited version 
- and derive new results related to the one-dimensional 
problem, in order to confront our approach to the avail- 
able exact results. In particular, we compute the corre- 
lation function, show that it takes a scaling form at long 
time and large distance and extract the associated scaling 
function. The NPRG results are in even better agree- 
ment with the exact results of Prahofer and Spohn [l| 
than those obtained under the MC or the self-consistent 
approximations . 

We stress the advantages of the NPRG formalism: 
(i) it is based on an exact fiow equation, which, given 
a microscopic model, yields the macroscopic properties 
of the system. The only input is the bare action, that 
is, no a priori knowledge, other than the microscopic 
model and its symmetries, is required to compute the 
physical observables. In particular, one does not have to 
assume scaling at long time and large distance nor the 
form of the scaling function; (ii) the approximations, 
which are mandatory to solve the exact fiow equation 
can be systematically improved and implemented in 
any dimension and, in contrast to most other nonper- 
turbative approaches, the calculation of higher orders 
is achievable and have been performed in practice in 
other systems [2^, [30] ; (hi) beyond scaling, a wealth of 
quantities can be calculated (for example, we compute 
here the correction to scaling exponent w); (iv) it has 
already yielded many nontrivial and accurate results 
in systems at equilibrium such as frustrated magnets 
31 1 , the random bond and random field _Ising model 
32 1, membranes [3j|, bosonic systems j34|, but also 
in systems out-of-equilibrium, where one can mention 
important advances in reaction-diffusion systems [35|. 

As the NPRG is a field theoretical method, our starting 
point is the field theory associated with Eq. ([1]), which 
follows from the standard procedure of Janssen-de Do- 
minicis relying on the introduction of a response field 
h and sources {j,j) [36]. This procedure allows one to 
explicitly carry out the integration over the Gaussian- 
distributed noise 77 in Eq. (|T|) upon doubling the number 
of fields (see e.g. [S^])- The generating functional reads: 

Z[j,j] - fv[h, ih] exp (^S[h, h]+ f jh + Th) (3) 



exponents in physical dimensions 12 1. 



S[h,h] = I \h[ dth-iyV'^h--{Vhf ) - Dh^ \{A) 



where x = {t,x). 

The remainder of the paper is organized as follows. In 
Section |lll we briefly review the nonperturbative renor- 
malization group formalism for out-of-equilibrium prob- 
lems. In Section Hill we analyze in detail the symmetries 
of the KPZ action ((4]) and derive Ward identities asso- 
ciated with the linearly realized ones. In Section |IV[ we 
build an approximation scheme based on a covarianti- 
zation procedure rooted in the symmetries, and derive 
explicitly an ansatz at the minimal order in the response 
field. The determination of the complete phase diagram 
and of critical exponents in all dimensions, using a sim- 
plified version of this ansatz, is reported in Section FVl 
Section |Vl] is then devoted to the calculation (using the 
full ansatz) of the scaling function and of some other uni- 
versal quantities in one dimension, which are compared 
with their exact counterparts. Technical details, such as 
some of the Ward identities, a discussion of the validity 
of our approximation scheme, the computation of vertex 
functions or the procedures for the numerical integration 
of the flow equations are reported in Appendices A, B, C 
and D respectively. 



II. THE NONPERTURBATIVE 
RENORMALIZATION GROUP 

The NPRG formalism relies on Wilson's RG idea, 
which consists of building a sequence of scale-dependent 
effective models such that fluctuations are smoothly aver- 
aged as the (momentum) scale k is lowered from the mi- 
croscopic scale A, where no fluctuations are yet included, 
to the macroscopic one k = where all fluctuations are 
summed over [23, [S^] ■ 

For out-of-equilibrium problems, one formally proceeds 
as in equilibrium, but with the presence of additional re- 
sponse fields and with special care required to deal with 
the consequences of Ito's discretization and with causal- 
ity issues, as stressed in detail in [S^l - from which con- 
ventions are taken throughout this paper. For future 
use, we define the Fourier conventions used in [38] and 
throughout this work: 



d'^p du 



(5) 



fit,x) 



(27r)'' (27r 
/(p) e-'P-^+"^* 



f{u,p) e-^P-^+-* (6) 
(7) 



where p = {io,p). 

To achieve the separation of fluctuation modes within 
the NPRG procedure, one adds to the original action 5, 
a momentum-dependent mass-like term: 



A5, = iy/i,(-q)K(q)],,/j,(q) 



(8) 



where the indices i, j G {1,2} label the fleld and response 
fleld respectively hi — h,h2 — h, and summation over re- 
peated indices is implicit. The matrix elements [Rf^{q)]ij 
are proportional to a cutoff function r{q^ / k^) (see Sec- 
tion [IIID]), with q = ||(7J|, which ensures the selection of 
fluctuation modes: r(x) is required to vanish as a; > 1 
such that the fluctuation modes hi [q > k) are unaffected 
by AiS„, and to be large when x < 1 such that the other 
modes {hi{q < k)) are essentially frozen. Since A5„ must 
preserve the symmetries of the model, we postpone the 
discussion of the precise structure of the matrix elements 
[Ri^{q)]ij to Section HlI Dl after the analysis of these sym- 
metries. 

In presence of the mass term AS^, , the generating func- 
tional ([3]) becomes scale dependent 



2n[j, 3hJ V[h, ih] exp ( -5 -AS^ + j jh + jhj (9) 

and the effective action r„[(p, (^], where ipi — {hi):-, are 
the expectation values of the fields hi in the presence of 
the external sources j and j, is given by the Legendre 
transform of Wk. = log Zk, (up to a term proportional to 

Rn) [Hii: 



T^[ip,Lp]+ log Z^[j,i] 



J.<y3j 



ipi [Rk\ij ipj- (10) 



From Fk, one can derive two-point correlation and re- 
sponse functions. 



[rK^^]ni2(xi,X2,(y5,(p) 



<5'r«[^,^] 



(5v3ii(xi)(5.^i2(x2) 



(11) 



and more generally n-point correlation functions that we 
write here in a 2 x 2 matrix form (omitting the depen- 
dence on the running scale k) 



.(n) 



(xi,...,X„,V3,^) = 



)"-^F(^)(xi,X2,y>,(^) 

S(pts{-K3)...S(pi^{-Kn) 



(12) 



The exact fiow for r„[(/3, 0\ is given by Wetterich's equa- 
tion which reads (in Fourier space) [29j: 



d^T^^^Ti I d^R^-G^ with G«: 



ri'^ + Rn 



(13) 



the full renormalized propagator of the theory. When k 
fiows from A to 0, F^ interpolates between the micro- 
scopic model Fk^a = S and the full effective action Ffj=o 
that encompasses all the macroscopic properties of the 
system [3^. Differentiating Eq. ((T5)) twice with respect 
to the fields and evaluating it in a uniform and station- 
ary field configuration (since the model is analyzed in its 
long time and large distance regime where it is transla- 
tionally invariant in space and time) one obtains the flow 
equation for the two-point functions: 

9.[F(2)],^.(p) = Try'9«i?(q).G(q).(^-ir(f(p,-p,q) 



.(3), 



.(3) 



+rf> (p, q) . G(p + q) • F ^^^ (-P, p + q) • G(q) (14) 



where the k and background field dependencies have 
been omitted, as well as the last argument of the F*^"^ 
which is determined by frequency and momentum 
conservation |38| . 

Solving Eq. ([T3l) (or Eq. ([T4|) ) is in principle equivalent 
to solving the model. In practice this resolution cannot 
be performed exactly since (J13p is a nonlinear integral 
partial differential functional equation. Hence one has to 
devise an approximation scheme. The main constraint on 
this approximation scheme is to preserve the symmetries 
of the problem. We thus now revisit the symmetries of 
the KPZ action. 



III. SYMMETRIES OF THE KPZ ACTION 

The KPZ action (|4]) possesses well-known symmetries, 
in addition to translation and rotation invariances: i) the 
Galilean symmetry and ii) the /i-shift symmetry, which 
can be expressed as the invariance of the action Q under 
the following transformations: 



h'(t,x) = X ■ V + h(t, X + Xvt) 
^' ^ h!{t,x) = h{t,x + Xvt). 

h'{t,x) = h{t,x) +c 



(15) 
(16) 



where v and c are arbitrary constant quantities. 

In one dimension, the KPZ equation also satisfies a 
fluctuation-dissipation theorem that fixes the exponents 
exactly. This property roots in a time-reversal symmetry 
of the action which, as shown in [39| , can be encoded in 
the transformation 



h'{t,x) = -h{-t,x) 

h'{t,x) ^ h{-t,x) + -^V^h{-t,x). 



(17) 



One can check that the invariance of the action (|3]) under 
the transformation (J17p requires that the contribution 
J V^/i(Vft,)'^ vanishes, which is 'incidentally' true only in 
d = 1. The time- reversal symmetry thus only holds in 
this dimension. This set of symmetries entails Ward iden- 
tities for the n-point vertex functions, as long as the mass 
term AS^ is appropriately chosen (see Section UlIPp . 

In fact, there exist even stronger symmetries of the 
KPZ action that, to our knowledge, were only pointed 
out in [40'|. They consist of 'gauging' in time the trans- 
formations i) and ii) in the following way: 



i') 



h'(t,x) ^ X ■ dtv{t) + h{t,x + Xv{t)) 
h'{t,x) = h{t,x + Xv{t)) 

n') h'{t,x) = h{t,x) + c{t) 



(18) 
(19) 



which we will refer to as Galilean-gauged and shift- 
gauged symmetries respectively. In these gauged ver- 
sions, v(t) and c{t) are arbitrary infinitesimal functions of 
time. Note that the action (j?]) is not strictly invariant un- 
der the transformations P^ and P^ but the correspond- 
ing variations of the action are linear in the fields, and 



this behavior also yields useful Ward identities, which 
we derive in the following. These stronger forms of the 
symmetries will be thoroughly exploited. 

Finally, an additional Z2 symmetry, which is manifest 
on the Cole-Hopf version of the theory, is nonlinearly re- 
alized in the KPZ action (j?]). The Cole-Hopf field trans- 
formation writes as follows: 



Mx) = ^logmx)| 
h{x) — w{x)w{x). 



(20) 



In terms of w and w - upon rescaling these fields and 
time - the KPZ action becomes: 

S[w, w]= j w{dtw - V^w) - -gb {wwf (21) 



with the bare coupling constant 

X^D 

9b ^ —^■ 



(22) 



This action is invariant under the simple Z2 transforma- 
tion (iv) w(t,x) — ?> iu{—t,x), ■w{t,x) — > w{—t,x). How- 
ever, in terms of the original fields h and h, this trans- 
formation becomes 

( h'it,x) = -hi-t,x) + -floglhi-t^x)] 
>\h'{t,x)^h{-t,x). ^^' 

Its highly nonlinear form renders complicated the in- 
duced Ward identities among vertex functions. Conse- 
quently, they are not given here as they will not be ex- 
ploited in the following. 



A. Shift-gauged symmetry 

The Ward identity associated with the shift-gauged 
symmetry can be derived by performing in the functional 
integral (jH]) the change of variables corresponding to the 
transformation (|19l) . As this operation must leave the 
value of the integral unchanged, one obtains that 

z{t)j{K) (h{^)^jdtc{t) {ASMt)Mm,j} 

(24) 
must vanish. Here, [_R„(q)]ii has been set to zero (as 
in Eq. (|4T|) ) for simplicity and for causality issues (see 
Section fHIPp . As the mass term is quadratic in fields, 
AS^lcit), /i(x)] is linear in h, and thus 

(A5.[c(t), h{^)])^j = ASMt), (H^)) - ^(x)]. 

The expression (jlOp of the Legendre transform then im- 
plies 



c(t)^^^-^(x)9*c(i)^=0. (25) 



After integrating by parts the second term, we conclude 
that the functional 



^Af^v] 



<^(x)at^(x) 



(26) 



is invariant under the transformation (1191) . In other 
words, the only non-invariant term J (pdtip of the bare 
action is not renormalized and the rest of the action is 
shift-gauged symmetric. 

Let us express this property on the n-point vertex func- 
tions, that we denote from now on as 

ri''")(xi,...,x,+„o 



which stands for the Tk~ vertex involving / (respec- 

tively m) legs - derivatives of F^ with respect to (p (re- 
spectively (p) - with the I first frequencies and momenta 
referring to the ip fields and the m last to the (p fields. 
At this stage, it is convenient to work in Fourier space. 
The 71-point vertex function are defined in Fourier space 
by 



(2^) 



d+l<:d+ 



i(^p,OF('-)(pi,...,p,+_o=/ r(' 



'(xi 



,x/+„)e 



J2iiSi-pi-U(^i) 



(27) 



xi---xi+„ 



r 



where again, the last frequency and momentum, fixed by 
translational invariance in time and space, are implicit. 
The shift-gauged symmetry entails that the n-point ver- 
tex functions in Fourier space satisfy the following prop- 
erty: 

F(™'")(l^1,P1 = 0, . . . , Pm+n-l) = ILOiS^lSnl (28) 

which means that, apart from the contribution of / (pdtip 
to Fk , the vertices vanish upon setting the momentum 
of one of the (p to zero. This is related to the fact that the 
field ip only appears in F^ with gradients. In particular, it 
roots the non-renormalization of the kinetic term (pdtp> of 
the KPZ action, which is well established in perturbation 
theory (see e.g. [l6|). 



B. Galilean symmetry 

1. Global Galilean symmetry 

Let us first review the Ward identities associated with 
the standard Galilean symmetry. As for the shift-gauged 
symmetry, one can prove that, as the transformation 
([T5|) encoding the Galilean symmetry of the bare ac- 



tion is affine, and provided the mass term AS^ is cho- 
sen Galilean-invariant (as in the form (J4T|) - see Section 
IIIIPp , F^ also possesses this symmetry. The correspond- 
ing Ward identity is: 



{x ■ v + Xtv ■ V(^(x)) 



Sp{:K.) 



+Xtv-Vp{-K)-— 



Sp{x) 



0, (29) 



that is Fk is invariant under the same Galilean transfor- 
mation (with the same parameter) as the bare action. 
One can then proceed to the derivation of the Ward 
identities for the vertex functions, by taking functional 
derivatives of Eq. (|29l) with respect to the fields and then 
evaluating them for uniform and static fields. For in- 
stance, one gets the standard identity for the three-point 
function: 



i— F^i)(c^ = 0,p = 0;c.i,pi) = Api,— Fii-i)(c^i,pi) 



dp 



duji 



(30) 

The Ward identity for a generic n-point function can be 
derived and one obtains: 



J 



*^4'"+i'")(w = 0,p = 0;pi;...;p,„+„_i) 



, , . 5 ^ d 

A ( Pl^ 1 \-p,n+7i-l^ 

OOJi OL0,n+n-l 



r 



F^")(pi;...;p™+„_i). (31) 



2. Galilean-gauged sym,metry 



the variation of the action under this transformation is 



Let us now come to the gauged form (fT8)) of the 
Galilean symmetry. As for the shift-gauged symmetry. 



linear in the fields and consequently, it entails a Ward 
identity which reads 



d''x\XVip{ic)^^~xdt ^^ 



6ip{x.) 



6ip{x.) 



6T, 



+AV^(x) ^^ - x9,V(x) 1=0. (32) 

From this functional identity, one can again deduce iden- 
tities for the vertex functions. They bare similar expres- 
sions as those for the global Galilean symmetry, but with 
a stronger content. For instance, the identity for the 
three-point function becomes: 

d 

iuj—rl^'^^uj^p:^ 0;uji,pi) = Xpi 

X (r^i)(^ + ^i,pi)-r^i)K,Pi)) , (33) 



which in the limit w ^ coincides with ([30| . The gauged 
identity ((33|) is stronger as it constrains the whole fre- 
quency dependence and not only the zero-frequency sec- 
tor. One could derive similar identities for generic n- 
point functions, but we now stress, instead, a more ef- 
ficient way to exploit the Galilean symmetry, that will 
guide our construction of the approximation scheme. 



3. Covariant time derivatives 

The previous Ward identities ensuing from the 
Galilean symmetry do not clearly refiect the geometrical 
interpretation of this symmetry. In order to do so, one 
can analyze the Galilean invariance from another angle. 
One can build quantities which are manifestly scalar un- 
der the Galilean transformation (jisp. upon introducing 
an adequate covariant time derivative. Let us define a 
function /(x) as a scalar under the Galilean symmetry if 
its infinitesimal transform under (1151) is given by 



<5/(x)=tAz7-V/(x). 



(34) 



With this definition, if / is a scalar then / df^xf is invari- 
ant under the Galilean transformation and can be used 
to build an action possessing the Galilean symmetry. It 
follows from this definition that the response field /i is a 
scalar, but that the field h is not unless one takes two 
successive space derivatives ViVj/i. The gradient of a 
scalar remains a scalar, but not its time derivative since 

5{dtf{^)) - \v ■ (tV(9i/(x)) + V/(x)) . (35) 

However, as in fluid mechanics, one can construct a co- 
variant time derivative 

bt = dt- AV/i(x) • V 

which conserves the scalar property, i.e. if / is a scalar 
then so is Dtf . Note that the covariant time derivative 



of the field h itself, that we denote -Dt, bares a special 
form with a 1/2 factor: 

Dth{^) = dth{^) - ^(V/i(x))2 (36) 

since h is not a scalar on its own. These covariant deriva- 
tives will constitute the building blocks in the construc- 
tion of our approximation scheme (see Section irvj) . 



C. Time-reversal symmetry in d = 1 

In d = 1, the action (|4]) with a mass term of the 
form (|^T|) exhibits the additional time-reversal symme- 
try. Note that this is a discrete symmetry (i.e. there 
is no infinitesimal transformation corresponding to (|17p). 
However, as it is linear in the fields, one can show using 
the same procedure as previously that F^; also possesses 
this symmetry. That is, it verifies: 

F,[</j(x),^(x)] - F4-</j(-i,x),(^(-i,f)+-^VV(-i,a;)]. 

(37) 
Again, one can derive Ward identities for the n-point 
vertex functions by taking derivatives of (1371) with respect 
to the fields and evaluating them at uniform and static 
field configurations. For the two-point functions, this 
yields, in Fourier space: 

F(^^")(..,p) = r(2-°)(-..,p) + ^p^T^y)(-^^p) 

+ ^p^Ti'.^)ico,-p) + [^yp^rrH~-,P), (38) 

which, given that F^ ' {uj,p) — F^ ' {—uj,p) and spatial 
parity reduces to: 



2ReF(M)(p)^_^p2r(o,2)(p), 

For the three-point functions, one gets: 

2ReF(i'2)(pi;P2) = --^p2r(o,3)(p^. p^). 
2ImFi2.i)(pi;p2) = -^pllmTi'-^\p,;p2) 



(39) 



D 



p?ImFli-2)(p2;pi). 



(40) 



-.(3,0) 



Similar expressions can be derived for F^ ' and for the 
four-point functions, which are reported for completeness 
in Appendix A. 



D. cutoff function 

The mass term AS^ defined by (|8]) must be chosen such 
that the Ward identities are preserved all along the flow. 
This is not obvious a priori. In particular, a quadratic 
mass term cannot be invariant under the shift-gauged 



symmetry. However, as explained above, if the varia- 
tion of AiSre under the transformation (|19p is linear in 
the field, the Ward identity remains identical even in the 
presence of such a term. Similar comments hold for the 
Galilean symmetry, which is preserved as long as the cut- 
off matrix R^, does not depend on the frequency. We thus 
choose a frequency independent matrix i?^ - which also 
preserves causality properties [38|. Moreover, the time- 
reversal symmetry imposes in one dimension the follow- 
ing relation between the components of the cutoff matrix: 
[Rii{q)\i2 = -ydI^ [-Rk(<z)]22, which we will impose in all 
dimensions. 

Finally, we set [RK,{q)\ii = since it is not necessary to 
consider a cutoff term proportional to (^(/j and since the 
presence of such a term would spoil the property that F^ 
is proportional to (^, important for causality issues (see 
[38| ) and also in order to preserve the shift symmetry. To 
summarize, the cutoff matrix is chosen with the following 
form: 



RM=r\^-, 







-2D, 



(41) 



where the running coefficients v^ and D^ are introduced 
for convenience [12| . Here we use r{x) — a/{exp{x) — 1) 
where a is a free parameter. 



IV. APPROXIMATION SCHEME 

Two routes have been followed in the past to build ap- 
proximation schemes to deal with exact NPRG flow equa- 
tions such as Eqs. (|13p and (|14p . Either one focuses on 
the long time and large distance properties, trying to de- 
scribe as accurately as possible the zero-momentum and 
zero-frequency domain while approximating the other 
sectors. Or the aim is to compute the momentum and fre- 
quency dependence of two-point functions and then the 
approximation concerns the three- and four-point func- 
tions. 

The usual implementation of the first route consists 
of performing a 'derivative expansion': F^ is expanded 
in powers of gradients and time derivatives of the fields 
|29l ISSf . This strategy has been undertaken to investi- 
gate the critical properties of many equilibrium [29] , and 
also noneq uilibrium systems such as reaction-diffusion 
processes [S^, where it has yielded satisfactory results. 
However, for the KPZ problem, the signature nonlinear 
term involves a gradient, and this seems to preclude the 
use of the derivative expansion [39| . 

The second route is the one followed here. At equilib- 
rium, one can compute the momentum dependence of the 



two-point function using the Blaizot-Mendez-Wschebor 
(BMW) approximation scheme which has proved to yield 
very accurate results for 0{N) models [41|. In the BMW 
framework, the momentum dependence of Fk and Fk 

(2) 

in the flow equation for F^ is truncated. However, a di- 
rect implementation of the BMW scheme for the KPZ 
problem is hindered by the symmetries which impose 
strict identities between F^ '" s with different m and 
at different momenta, see e.g. Eq. pip , and render com- 
plicated the writing of a truncation. 

To overcome this difficulty, we propose, instead of de- 
signing directly truncations on the three- and four- point 
functions, to construct an ansatz for F^ that manifestly 
preserves the shift-gauged and Galilean-gauged symme- 
tries while allowing for arbitrary frequency and momen- 
tum dependencies in the two-point functions. This is 
easily achieved by combining together the basic Galilean 
scalars of the theory - ip, "S/iS/ jip and Dtf, and arbitrary 
powers of their gradients and covariant time derivatives 
Dt - since, by construction, all calculations involving a 
functional of these quantities automatically satisfy the 
Galilean symmetry (the other symmetries, but the non- 
linearly realized one, being easily enforced a posteriori). 

On this basis, various truncations may be performed 
to obtain an ansatz for F^. One could, for example, 
keep complete dependencies in (p and V^(^, while treat- 
ing only the zero-frequency sector. This 'super derivative 
expansion' option is left for future studies. An alterna- 
tive choice, pursued here with the ultimate aim of con- 
fronting the obtained two-point correlation functions in 
one dimension with the exact ones, is to truncate the field 
dependence of F^ while preserving the complete depen- 
dence in both momentum and frequency of the two-point 
functions. 

In the following, we hence proceed to a field truncation 
at the minimal order. It appears that the truncations of 
F^ at a given order in either ip or tp are not equivalent 
because the dependence in p, contrary to that in ip, is 
constrained by the Galilean symmetry - p enters in Dt. 
The minimal order in fields amounts to keeping terms, (i) 
at most quadratic in ip, (ii) linear in V^(^ or in Dtp, and 
(iii) to combine them with arbitrary powers of Laplacian 
V^ and covariant time derivatives Dt. This choice allows 
for an arbitrary dependence in momentum and frequency 
of the two-point functions. Notice that the dependence 
in p is not restricted to be polynomial since arbitrary 
powers of this field are included through the covariant 
derivative Dt. To summarize, the ansatz for F^ consid- 
ered in the following reads: 



J 



r^[p,p] = 



pf^{-Dl -V')Dtp - pf^{-Dl -V')p - — [v\rA-Dl -V')p + prA-Dl -V')V'p 



(42) 



Setting /^ = 1 and f^ ~ f!^ ^ D in the previous ansatz, 
one is left with the bare action (j4]) - written in terms 
of covariant time derivatives. Allowing non-constant 
(and scale-dependent) functions f^ {X = D, v, A) of the 
Laplacian and of covariant derivatives is then the way to 
include an arbitrary dependence in momentum and fre- 
quency for the two-point functions that preserves all the 
symmetries, at the lowest (second) order in ip. 

The operator Dt is squared in /^ since this function 
is real. For /^ and /^, to take the square of Dt is a 
choice that ensures a natural separation between real and 
imaginary parts of Fk ' , without loss of generality for 
the two-point functions. Note that one needs to make 
precise the meaning of the expressions /^(— iPf,— V^) 
{X = D, V, A). Indeed, it is ambiguous as it stands since 
the Dt and V operators do not commute. We take the 
convention that all Dt 's are on the left of the Vs. We also 
assume that these functions can be expanded in series of 
their arguments: 



f^i^bi-y^)^ J2 a,nn{~Dir{~y^r 



(43) 



m,n— 



One can then consider the constraints stemming from 
the other symmetries. The shift-gauged symmetry im- 
poses that /^(w^,p^ = 0) = 1. In dimension one, the 
time-reversal symmetry implies that the two functions 
f^ and /^ become identical, and also that f^{uj^,p^) = 
1. In a generic dimension d, the ansatz ((42]) thus consists 
of three independent running functions oip^ and w^, and 
is reduced to a unique function of p^ and w^ in d = 1. 

To compute the flow equations of the functions 
/^(— i)^, — V^), X = X,^ 01 D, one needs the expres- 
sions of the n-point functions up to n = 4. All the 
calculations are performed at vanishing fields. For the 
two-point functions, they are straightforward and yield 










:f^{u;^p^) + ^p'f;:iu:^p')iU) 



vt:{^\r). 



In this case, there is no ambiguity arising from the or- 
dering of the gradient and covariant time derivative op- 
erators. The calculation and the explicit expressions of 
the three- and four-point functions, more lengthy, are de- 
tailed in Appendix C. One can check that all the Ward 
identities derived in Section IIIII are satisfied by these 
functions. 

We can be a priori confident in our approximation 
scheme, as it is close in spirit to the BMW approximation 
which can accurately capture the momentum dependence 
of two-point functions 41]. A more detailed discussion 
on the validity of our approximation scheme can be found 
in Appendix B. 



V. SIMPLIFIED APPROXIMATION 

The ansatz ((321) in generic dimensions d can be fur- 
ther simplified by restraining the form of the functions 
/^(— I)j, — V^). The idea is to focus on a reliable de- 
scription of the zero-frequency and zero-momentum sec- 
tor of the theory while circumventing the limitations of 
the derivative expansion which is problematic here. For 
this, one can neglect all D^ dependence in these func- 
tions: f^i-Dl-V^) -^ /£i--V^) for X = A, z/ and D. 

The ansatz proposed in [12l | - which was derived be- 
fore the formalization of the Galilean-gauged symme- 
try in terms of covariant derivatives - is actually en- 
dowed with an additional simplification since we im- 
posed the one-dimensional constraint /^ = 1 in all di- 
mensions (or stated otherwise we extended the condition 
f^{p — 0) = 1 holding for all d to the whole p sector). It 
reads 



T4^,^] = / {^(x)A¥'(x)-^(x)/,^(-V2)^(x) 

- ;^^(x)/:(-V2)vV(x)} . (45) 

As a matter of fact, this ansatz implies that all interac- 
tion vertices (F^ ■^ with n > 2) are reduced to their bare 

(2 1) 

form. Among these, the only non- vanishing one is F^ : 
^';^'^H^i,Pi,^2,P2) = Xpi-P2- (46) 

With the ansatz ((45)) . the two-point functions in generic 
dimensions simplify to 



r(^'"Hc^,P)^o 
ri'-'H^,p)-^o. 

Fp)(c.,p)=-2/f(p2). 



F(l■l)(..,p)=^.. + -^pV:(p2) (47) 



With the aim of analyzing the fixed point structure, 
one introduces dimensionless and rcnormalizcd functions 
/k = f^/'^n (respectively /f = f^/D^), where i/„ (re- 
spectively Dk) are running coefficients (identifying at the 
microscopic scale k — A with the bare parameters i' 
(respectively D) of the action ((H)) that are related to 
anomalous dimensions Kd^i^K — — '7iy('t)'^K (respectively 
ndf^Df^ = —ri,:,{K)Di^). At a fixed point, these coefficients 
are expected to behave as power laws v^ ^ n^^" and 
Z?K ~ K~^" and the scaling exponents are then expressed 
in terms of these anomalous dimensions as 2; = 2 — ry* and 

The flow equations for the functions /^ and /^ can 
be found in [12|. The Galilean invariance ensures that 
the flow of the dimensionless running coupling constant 
g^ = k'^~^ X^ D K, / ly^ is reduced to its dimensional part 



ds9K ^ gnid- 2 + 377^(k) ~ rio{K)) 



(48) 



with ds — kOk , which thus enforces the identity z + x = 2 
at any non-Gaussian fixed point. The subleading expo- 
nent uj, governing the corrections to scaling, is indepen- 
dent of the leading critical exponents x and z and can be 
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FIG. 1: (Color online) (a) Flow diagram stemming from our 
(simplified) minimal order approximation in the (d, gvd) plane 
(with v^^ = 2'*+^7r''/^r[|] a normalization factor related to 
the integration volume). Red circles: renormalized value g%,Q 
at Fsc. Dashed purple line: Gaussian fixed point Few- Cyan 
squares: bare value g^, which separates the basins of attrac- 
tion of Fsc and Few. Gray arrows symbolize flow lines, (b) 
Variation with d of x = 2 — z for fsc (red circles: our re- 
sults; orange squares: numerical values from [43l . |44J . See 
also Table [2) 



mates for the roughness exponent x and the subleading 
exponent w. For d ^ ^^ x decreases almost linearly with 
d, with the exact value \ (respectively |) recovered in 
d=\ (respectively d = 0), and a reasonable but deteri- 
orating agreement with numerical values for 2 < d < 4. 
In higher dimensions, x increases with d, at odds with 
both numerical values and the scenario of d = 4 being an 
upper critical dimension beyond which x = [19, 20, 45]. 
The values for the critical exponents can be refined by 
computing the next order of our approximation scheme, 
which is work in progress. 

Regarding Frt, we record negative values of x for 
2 < d < 5, which is reminiscent of perturbative results 
performed at fixed d |14| but in contradiction with exact 
results dictating x = [lE, |4a] • This is to be imputed to 
our breaking of the Z^ symmetry manifest in the Cole- 
Hopf formulation and non-linearly realized in the KPZ 
problem. 



VI. ONE-DIMENSIONAL SCALING FUNCTION 



TABLE I: Exponent values in integer dimensions. The aver- 
age numerical values are extracted from 43]. To our knowl- 
edge, no estimates of u) are available in the literature. 



d 


1 


2 


3 


4 


X [38] 


0.50 


0.33 


0.17 


0.075 


X (numerics) 


0.50 


0.38 


0.30 


0.24 


uj [38] 


0.817 


0.70 


0.63 


0.54 



calculated from the flow behavior as s — >• — oo, using for 
instance that 



We now consider the full ansatz P^ to calculate the 
momentum- and frequency-dependent two-point correla- 
tion function of the one-dimensional problem and to ex- 
tract from it the universal scaling function. We recall 
that, due to the incidental time-reversal symmetry, the 
ansatz (|42p simplifies in d = 1 to only one running func- 
tion left: f^ = f^ = /k and /^ = 1. We also drop from 
now on the vector arrows and set u = D = 1 since these 
two coefficients can be absorbed in the action ^ through 
an appropriate rescaling of the fields and of time and the 
change of coupling constant A — > y/gl — XD^I"^ jv^^"^ . 



Ok^ g* + 51 exp(ws) 



(49) 



in this limit. 

Despite its simplicity, this ansatz yields the correct 
phase diagram in all dimensions (Fig. [Ija)) - including 
the strong coupling fixed point, and reasonable exponent 
values in physical dimensions (Fig.[ljb)) [12]. More pre- 
cisely, in all dimensions studied (i.e. up to d = 8), there 
exists besides the Edwards- Wilkinson fixed point Few, 
a fully- attractive non-trivial strong-coupling fixed point 
Fsc- In all d, generic scaling is found, i.e. the flow al- 
ways reaches one of these fixed points. For d < 2, Fgc is 
reached from any initial condition. For d > 2, Few f-nd 
Fsc become locally fully attractive. More precisely, there 
exists a critical value of the (bare) coupling, denoted gc, 
that separates the basins of attraction of Fsc and of Few 
and for which the flow reaches the roughening transition 
fixed point Frt- In d = 2, Frt coincides with Few, a-nd 
becomes non-Gaussian for larger dimensions. 

Because of Eq. ([48]) and since the fixed point coupling 
constant gg^ is non zero, Fgc is characterized by a single 
exponent. We thus only discuss below the values ob- 
tained for X UM ■ Table U and Fig. [IJb) contain our esti- 



A. Flow equations 



The flow equation for the function /„ can be obtained 
either from the flow of F^ ' ' or from the one of T^ ' 
according to (H^ . These two flows are equal in d = 1 
since the ansatz preserves all the symmetries including 
the time-reversal one. They can be computed using Eq. 
([T4|) . where Gk is the propagator matrix G^ = [Tk + 
Rk.]~^ as defined in Eq. (fT3|l . Using the expressions (|4T|) 
and ([m (specialized to d = 1, i.e. setting i/^^ = £)„, 
Ik — Ik = Ik, ^ud /^ = 1) for the matrix elements of 

.(2) 



i?K and Fk ' , one finds for the propagator 



GK(w,g) 



2fc,(L.2,q2) Y,iL,,qy 



P.iu^,q^) \ Y*{u,q) 







(50) 



where K{uJ^,q^) = f^{u:'',q^)+D^r{q^/n'), Y^{LO,q) - 

iuj + q^k^{uj^,q^) and P^{uj^,q^) = lo^ + {q^k^f . The 

calculation of the flow equation for e.g. F^ , involving 
matrix products and trace according to Eq. (|14p . is then 
straightforward. This yields for the flow of /«: 
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dj,iw,p) = -\j 0«5«(<z)|-p2(^i-^X«(c.2,g2)r(^2,2)(q^_q^p) 



X 9 



y«(tn + 0.,^ + g)r(,i'2)(_q, _p) + 2fc,((tn + u)\ {p + g)2)r(,24)(p + q, _q)^ 
(y;(tn + c.,p + q)Ti'''^ (q, p) + 2fcJ(tn + ojf, {p + qf)r^^''^ (q, -p - q)) 
^(r;f(a.,g)r(i'2)(p + q,_p)_2Xj..2^q2)r(2,i)(p + q^_q) 



(51) 



where q = (w,q), 



P = (^,P), ^K(w^g2) = w^ - 
q^kf^{u)'^,q^), S^{q) = Df^r[q^ /k^) and where the expres- 
sions for the three- and four-point vertex functions Y^ 
are given in Appendix C. The flow equation for /„ is 
hence an integral over the internal momentum q and fre- 
quency u and depends on the external momentum p and 
frequency w. 

As we intend to study the fixed point properties, we in- 
troduce dimensionless and renormalized quantities. Mo- 
menta are measured in units of k, e.g. p — p/k, and fre- 
quencies in units of Z^^k^, e.g. U7 — zu/^Df^n'^). We define 
the dimensionless renormalized function /« = f^/D^, as 
at the bare level, /k=a = D = 1. The sole running 
anomalous dimension 77^ is defined by Kd^^lnD^ — —^n 
so that Z?K ~ n~^ at the fixed point. The critical expo- 
nents in d = 1 are then given hy z = 2 — if and X = r]* . 
The absolute normalization of /« and D^, is fixed by set- 
ting /k(0,0) — 1 for simplicity ^4]|. 

The flow equation for the dimensionless function reads: 



B. Numerical integration and error bars 



9s/k(^,p) = d, 



~ Jk. 



1 

ri^f^{w,p) 



D^K^ 



(2 -rin)'!^ d^fK{w,p) 
p dpfjw,p) + jr-Kd^fi^{m,p) (52) 



where dg = nd^,. Once the substitutions for dimension- 



less quantities have been performed in ((5T|) . the last term 
■^dsJK. in (EH) is dimensionless, and depends on the ex- 
ternal dimensionless momentum p and frequency w. 

The dimensionless running coupling constant reads in 
one dimension ^„ = k~^ X^ / D^. Its flow equation is again 
reduced to its dimensional part due to Galilean invari- 
ance: 



ds9K = gK(2?7K - 1) 



(53) 



and one finds as expected that X — V* — ^ ^^ any fixed 
point with g* ^ 0. The subleading exponent uj is calcu- 
lated according to (^^ and we find values in agreement 
with the ones obtained with the simplified ansatz of Sec- 
tion |V] (see Tables |T] and HH) . 



The numerical integration of Eqs. (|5^ and ([55)1 was 
performed using standard techniques and the numerical 
error stemming from the integration procedure was as- 
sessed by resorting to several system sizes, resolutions, 
etc. (see Appendix D). It turned out that the numeri- 
cal error was negligible compared with the error coming 
from the approximation scheme itself, which constitutes 
the dominant source of inaccuracy and which we now 
discuss. 

Ultimately, the error relative to the whole approxima- 
tion scheme can be evaluated by comparing successive or- 
ders of approximation. However, when only a given order 
is available, one can assess the accuracy of this approx- 
imation by varying the parameter a of the cutoff func- 
tion. Indeed, although physical quantities are obtained 
in the limit k — > (where Rk=o vanishes) from the exact 
NPRG fiow equations ^^ or ([H]) and thus do not de- 
pend in principles on the choice of the particular profile 
of i?K, any approximation introduces a spurious residual 
dependence on the regulator [3(1 . This sensitivity to the 
cutoff can be investigated and exploited to estimate the 
accuracy of the approximation implemented. 

Here, a was varied between 0.5 and 14, and all the 
physical quantities computed and discussed in Section 
IVICI exhibited a similar large plateau behavior. The 
value given in the following (text and Table |TT| for a 
quantity corresponds to the central value on the plateau, 
and the error bars reflect the dispersion of this quan- 
tity around its plateau value when a varies in the range 
[0.5, 14]. We emphasize that this procedure only provides 
an intrinsic error at a given order of approximation re- 
lated to an artificial regulator dependence. Accumulated 
experience in NPRG calculations seem to indicate that 
in general this procedure yields a reasonable estimate of 
the order of magnitude of the committed error [30j . How- 
ever, in some cases, these error bars may represent only 
lower bounds on the uncertainties. Definite error bars 
can be inferred from the comparison with the next order 
of approximation. 
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FIG. 2: (Color online) Typical shape of the dimensionless 
fixed point function f*{-&,p) (recorded at s = —25 and for 
a = 1). 

Runs were started from various different initial condi- 
tions including the bare action (/a = 1) with ^a typically 
between 1 and 10. In all cases, the function f,^ smoothly 
deformed from its initial shape to eventually reach a fixed 
point where it stopped evolving, typically after s < —20. 
In other words, the NPRG flow equations entail generic 
scaling for the one-dimensional KPZ problem, since the 
long time and large distance behavior is always governed 
by this non-trivial fixed point, yielding a scaling regime 
(described in Section IVICp . Fig. [2] shows a typical shape 
of the fixed point function f*{m,p). It is fixed by the 
normalization condition to unity at vanishing momen- 
tum and frequency and decays as a power law when ti7 
and/or p become large. 



C. Scaling function 

1. Extraction of the scaling function 

We now turn to the description of the fully attractive 
fixed point. At this fixed point dsfn{'&,'p) = (in Eq. 
([521) ') and g* 7^ such that rJK = rj* = \. Moreover, 
we verified numerically that the nonlinear term dsfi^/D,^ 
in ([52]) decouples (that is dsj^l^K — > 0) when w ^ 1 
and/or p 3> 1. As a consequence, at the fixed point, 
Eq. (152]) reduces in the regime w and/or p ^ 1 to the 
homogeneous equation 

\r(w,v)+vdir{^.V) + \^ d^f*{^,p) = 0. (54) 
Its general solution has the form 



r{^,p) 



1 



1/2 



p 



w 



3/2 



(55) 



where the function C, cannot be determined from the ho- 
mogeneous equation but can be extracted from the nu- 



merical solution of the full equation ((52)) - by tabulating 
the values p^^'^f*{w,p) against the ratios vjjjp'l'^. We 
observe that the fixed point function /*(ro,p) is regular 
for all vj and p (see Fig. [2]), which allows us to deduce 
the limits of C- First C(0) has to be finite for the limit 
TO —J' to exist and followingly /* ^ C(0)P~^^^ for p > 1 
at fixed 717. Second, C must behave as ^(2;) ^ CooS^"^ 
as X — )■ 00 for the limit p — >■ to exist and followingly 
/* ^ Coo^"^^'^ for K7 :^ 1 at fixed p. We show below 
that the form of the solution ([55)1 entails scaling for the 
(dimensionful) correlation function C(n7,p). 
We now consider the fixed point dimensionful function 



fim,p)^D^f*{m,p)^DJ* 



P 



D^K^ 



(56) 



The physical limit is obtained when the running scale 
K tends to zero at fixed values of w and p. This limit 
is precisely equivalent to to » 1 and/or p ^ 1, which 
corresponds to the scaling regime described above where 
/* takes the form ([55]) . Moreover, when k — > 0, D^ 
behaves as a power law D^, ~ Dqk^^^^ where Dq is a 
non-universal constant. Hence, the physical dimensionful 
function / is expressed in terms of the scaling function C 
as 



/(^'P) = ^^ 



1 



Dqp 



,3/2 



(57) 



The correlation function is defined, upon inverting the 
matrix F*^^' of the two-point vertex functions, by 



r(°-^)(TO,p) _ 2/(to,p) 
|r(i.i)(TO,p)|2 w^ + p*fiw,py 



(58) 



Jk — /' Jk — ^' 



C(to,p) 



where the ansatz (|44)l . with /^ 
and 1^ = = 1, has been used in the second equality. 
Note that the definition ([55]) of the correlation function 
corresponds in real space to the connected mean value 
{h{t,x)h{0,0))c which differs by a factor —2 from the 
definition ^ in the co-moving frame. Restoring the v 
and D coefficients also yields an additional factor, such 
that the correlation function defined by ([58]) finally re- 
lates to the Fourier transform of the correlation function 
defined by ^ via an overall multiplicative factor 



Co 



2 D' 



(59) 



Replacing in (|58p the dimensionful fixed point function 
/ by its expression ([57]) yields 



C{m,p) 



D, 



0'' \Do p3/2 



57/2 



r^yp^ + D^^^^^ 



^'/' Do ^2/^3 + ^"2 (_0^) 



1 



-F 



pV^Do V^op3/2 



(60) 

(61) 
(62) 
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FIG. 3: (Color online) Leftmost curve: Scaling function 
^[^W^l corresponding to the data collapse ((62}, with the 
two insets zooming different parts to show the very high qual- 
ity of the collapse (note the vertical scales); the various colors 
(gray shades) differentiate the contributions of distinct values 

of ro. Rightmost curve: Scaling function F ir = -:Sj 1 mul- 
tiplied by t'^'^ to illustrate the power-law behavior of the tail 
F ~ T-~'^/3 (see text). Note that the small dispersion visible 
in the inset zooming the tail for r > 100 is hence magnified 
by a factor of order 100^' "^ . 



Eq. (1501) shows that the correlation function takes a scal- 
ing form at large distance and long time. Let us em- 
phasize once more that we did not assume scaling, the 
latter emerges from the fixed point solution of the fiow 
equation starting from any reasonable microscopic initial 
condition. The scaling function F is hence universal, up 
to an absolute normalization and a rescaling of its argu- 
ment by a multiplicative constant. 

According to Eqs. ((6T|) and (|62|) . this function F can be 
computed from the numerical data for the dimensionless 
fixed function, selecting its values for arguments in the 
range w and/or p ^ 1 which corresponds to the scaling 
regime. We proceeded for various initial conditions and 
for different values of the parameter a. In all cases, we 
observed for each value of a the expected data collapse 
- generating the scaling function F - with a very high 
precision, as illustrated in Fig. [3l 



2. Normalizations 

Our aim is now to confront our scaling function F with 
the exact result obtained in Ref. pj. The first step con- 
sists of fixing the normalizations. Indeed, two constants 
are involved in the standard definition (as recalled be- 
low Eq. ^ ) of the universal scaling function g from the 
correlation function C 

C{t,L) = const.T2x/^g(const.'L2xT-2x/^)^ 



that have to be fixed. The scaling function g{y) is nor- 
malized in jll in the following way: 



9iy) 



lim 



c({2XAf 



,1/3 



y,t 



[IXAH) 



2/3 



(63) 



with A = D/v. Then three functions /, / and / are 
introduced in [l|, through the definitions 



I{y) = -^9"{y) 



f{k) = 2 dy cos{ky)f{y) 
Jo 

/•OO 

/(r) = 2 / dk cos(fcT)/(fc2/3) 



(64) 



and imposing the additional normalization condition 

/(O) = 1. (65) 

The last function / is proportional to C(a;,p)p^/^ (see 
[i|), that is it corresponds to our function F recon- 
structed using (|62|) up to normalization factors. The 
analog of the functions / and /, that are computed in 
the following, will also be denoted with capital letters F 
and F respectively. The precise normalization constants 
between F and / can be established as follows. First one 
deduces from Eqs. (|63p and ([M)l that: 



fir) 



n7/2 



rC 



n3/2 



25/3_X4/3^5/3i7/3 \l^2XAt^)^^^ ' * 



(66) 

Then, comparing this expression with our definition (j62p 
and taking into account the multiplicative factor (1591) , we 
obtain the relation 



fir) 




(67) 



where we used that the bare value A is related to the fixed 
point coupling g* through g* = gb/D'^ = ^^D/{DqI^^). 
Relating the Fourier transforms is then straightforward: 




(68) 



m-2^F 



Finally, using the same normalization criterion (j65l) as 
in [l| gives the absolute vertical normalization _F„ = 
j^ = F{0) and the absolute horizontal normalization 



(4F„5*) 



1/3 



3. Properties of the scaling function f 

Let us first compare the scaling functions / and F. 
The function / is studied in detail in [l| and also in 



[19|,l23|, where interesting features are highlighted. Ac- 
cording to Ref. [l[ , the function / first decreases to vanish 
at ko ~ 4.36236 . . . and then exhibits a negative dip of 
coordinates (fcrf ~ 4.79079 . . . Jd ^ -0.0012023 . . . ). Af- 
ter this dip, the function decays to zero with a stretched 
exponential tail, over which are superimposed tiny oscil- 
lations around zero, only apparent on a logarithmic scale. 
A heuristic fit of this behavior for /c > 15 is given in [l| 



f{k) - 10.9fc- 



-9/4 



sm 



^3/2 



1.937 e~2 



1,3/2 



(69) 



We show below that we here recover qualitatively all 
these features, with reasonable estimates for the different 
parameters that characterize them. Following (|64)) . the 
function F is defined from F by the integral 



F{k) 



dr 



-i{Tk^/'^)F{i 



(70) 



which has to be computed numerically. The function F 
stems from the superposition of the numerous curves in- 
volved in the data collapse of Eq. (p^ . Although the 
data collapse is excellent, there subsists a high frequency 
and small amplitude noise due to the re-ordering of the 
points (apparent with a large zoom as presented in the 
insets of Fig. [3]) . As we intend to investigate the proper- 
ties of the tail of F with high precision, we first need to 
eliminate this noise. For that, we devise an appropriate 
family of analytical fitting functions to smooth our data 
for F{t). The choice of the family of fitting functions is 
determined as follows. 

First, they have to reproduce the large r behavior of 
F. The latter can be inferred from the limits of the 
function C,{x) established after Eq. (I55|) . One obtains 
F{t) ^ FooT~'^/^ as T — ?► oo. For each value of the param- 
eter a, the proportionality constant Fqo can be estimated 
graphically (up to 4-5 digits) from the curve t'^/'^F{t). 
Moreover, the fitting functions have to be even, and fi- 
nite at the origin. To satisfy these constraints, we build a 
family of fitting functions as an expansion in elementary 
rational polynomials of r^ raised to the adequate power 
to reproduce the large r behavior: 



Ffit(T) 



Qqo + oo2r^ 
1 + aoiT^ + ao3T** 



Oio 



1- 

020 



flllT^ 



1 + a2\T^ 



7/6 



(71) 



with the additional constraint (002/003 -I- oio/on -I- 
020/021 + ...y'^ = Foo m. The first fit is achieved 
using the three independent coefficients a^i. The corre- 
sponding denominator introduces four non-analyticities 
in the points of the complex plane ±Zq with 



2a, 



03 



-ooi ±i\ -a, 



^01 
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(72) 
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FIG. 4: (Color online) Illustration of the fitting procedure 
described in the text, using the same data for F as in Fig. 
|3] but with unified colors (black curves: on the left F{t); 
on the right F{t)t'^'''^). The four successive orders of the 
fit appear superimposed in the main graph; the two insets 
present large zooms (note the vertical scales) to differentiate 
them and highlight the convergence: the higher the order, 
the closer to the data (black-dotted curve) the fit lies. Note 
that for the right curve and inset F{t)t ' the differences are 
again magnified by a factor r^". 



whose coordinates were found to be very robust against 
the choice and the order of the fit as explained below. 
Note that these four points appeared never to lie on the 
imaginary axis, nor on the real axis (— Oqi + 4ao3 > 0). 
We then add in turn simple monomials corresponding 
to purely imaginary poles up to nine independent coef- 
ficients (we observed that additional complex poles were 
systematically decomposed into two purely imaginary 
poles). This procedure turned out to converge rapidly 
for all a values, as illustrated in Fig. |4l In particular, 
the existence and the coordinates of the complex singu- 
larities ±z^ appeared to be a robust feature of the data. 
Moreover, they are always found to lie the closest to the 
real axis when other poles are included according to (j7ip 
|48| . For the remaining calculations we work with the 
analytical expression Ff^t which reproduces faithfully our 
data, and hence drop in the following the index gt. 

We computed numerically the function F from F ac- 
cording to (|70p and normalized it using (|68p . which 
yielded the following results. We find that F exhibits 
all the qualitative features of the exact function /, in 
particular the existence of the negative dip and the 
subsequent stretched exponential decay with the pres- 
ence of oscillations in the tail, as displayed in Fig. [5l 
Moreover, it compares very accurately with the exact 
function. Indeed, we find for the position of the first 
zero fco — 4.21(4) and for the coordinates of the dip 
{kd — 5.07(6), Fd ~ —0.013(1)), which are reasonably 
close to the exact results, see Table HIl As for the behav- 
ior of the tail of F, it can be inferred analytically from 
the pole structure of F. Let us define F{k) = F{k^^^), 



which is hence the standard Fourier transform of F 

F(fc) = / — cos(Tfc)F(r). 
Jo '^ 



(73) 



As F{t) is C°° , the tail of F is dominated by the sin- 
gularities of F{t) in the upper complex half plane lying 
the closest to the real axis, which are Zq and —Zq . De- 
noting Zq 



ingly 



_j_ gikz, 



= ao + ibo with bo > 0, we obtain F{k) ^ 
ex e"^"*^ cos(aofc) as fc -> oo and foUow- 



F{k) ^ F{k^/^) 



e-""''"'" cos(aofc3/2) 



as fc — ;> CX3. 



(74) 

F hence decays following a stretched exponential with 
superimposed oscillations on the scale k^^^ exactly as ob- 
served in the exact solution (see ([5^). Regarding the tiny 
magnitude over which develop these features, this agree- 
ment is remarkable. Let us emphasize that within the 
MC approximation and the self-consistent expansion, the 
stretched exponential behavior and the oscillations were 
also highlighted, but not on the correct scale [l|, [l^ [2j] . 
The NPRG method hence seems to provide more accu- 
rate results than other nonperturbative approaches, and 
with much less inputs (no assumed scaling nor a precise 
scaling form). 

A consequence is that, conversely to the MC approach 
for instance, the coefficient of the exponential and the pe- 
riod of the oscillations can here be estimated - extracted 
analytically from the values of zq- We find bo — 0.34(1) 
for the coefficient of the exponential, and Oq — 0.17(1) for 
the pulsation of the oscillations (checking that these val- 
ues coincide with the same quantities estimated graphi- 
cally from the tail of F computed numerically - the latter 
being determined with much less accuracy). These coeffi- 
cients can be compared with the coefficients of the fit (p^ 
of the tail of the exact function which are 6o = Qo = o- 
They are of the same order, which is already highly non- 
trivial. The discrepancy is visible on a logarithmic scale, 
as illustrated in the inset of Fig. [SJ Let us emphasize once 
more that the error bars attributed here to a given quan- 
tity reflect the weak variations of this quantity around 
plateau values when the a parameter of the cutoff func- 
tion is varied. It represents an error estimate intrinsic 
to the order of approximation under study and does not 
imply that the value at the next order of approximation 
would necessarily fall within these error bars. 



4-. Scaling function in real space and universal amplitude 
ratio 

We finally come to the real space scaling function / 
defined in (j64p . We computed it numerically by Fourier 
transforming F 



1'°° dk 
F{y) = / — cos{ky)F{k). 




FIG. 5: (Color online) Comparison of the scaling function 
F{k) (red curve with dots) obtained in this work with the 
exact one f(k) (black curve with squares) from [l|. The inset 
shows the stretched exponential behavior of the tail with the 
superimposed oscillations, developing on the same scale fc^". 
Note the vertical scale: this behavior develops with ampli- 
tudes below typically 10~ (see text for a detailed comment 
of the figure). 



It is compared in Fig. [6] with the exact result. Again, the 
function F reproduces accurately the exact one /, though 
it is the less accurate of our three scaling functions since 
it stems from the two successive numerical (oscillating) 
integrations (|70|) and ((75|) of our raw data. The tail is 
particularly sensitive to this loss of precision. The exact 
function /(r 



(y) is 
0,^ 



IS found to decrease as exp(— cy ) when 
2/ — )■ 00 m [ijj, whereas the decay of the function F{y), 
though it starts with the correct behavior, rapidly crosses 
over to a simple exponential decay exp(— c'y). 

From this function can be computed the universal am- 
plitude ratio go with the definition given in [l| 



50 



dyyfiv)- 



(76) 



(75) 



The exact result is the Baik-Rain constant 3 50 = 
1.1503944783... (more digits can be found in [1|). Here, 
the universal constant g^ can be estimated by perform- 
ing the numerical integration corresponding to (|76|) using 
F . However, this amounts to achieving three successive 
numerical integrations from our raw function F and the 
resulting precision is low. Alternatively, one can compute 
part of the involved integrals analytically, using the def- 
inition given in [ij: go = 5(0) where g{y) is the original 
scaling function proportional to the second derivative of 
J{y) (see ([M]) ). Indeed, g{y) can be expressed integrating 
twice / through its Fourier transform / 

g{y) — .g(0) = — dv du dk cos{ku) f(k) 



(77) 



^1"* (l-cos(fcj/))4?- 
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exact 



NPRG 



FIG. 6: (Color online) Comparison of the scaling function 
F{y) (red curve with dots) obtained in this work with the 
exact one f{y) (black curve with squares) from [l|. 



go 


1.15039 


1.149(18) 


feo 


4.36236 


4.21(4) 


kd 


4.79079 


5.07(6) 


u 


-0.00120 


-0.013(1) 


ao 


1 

2 


0.17(1) 


bo 


1 
2 


0.34(1) 


LJ 


- 


0.82(4) 



TABLE II: Characteristic parameters of the different scaling 
functions, from the exact results of Ref. [t] and from this 
work: i) relative to /, universal amplitude ratio go; ii) rel- 
ative to /: position of the first zero fco, coordinates of the 
negative dip {kd,fd), coefficient of the stretched exponential 
&o, pulsation of the oscillations ao; iii) correction to scaling 
exponent uj. The error bars reflect the weak variations around 
plateau values when the Q-parameter of the cutoff function is 
varied between 0.5 and 14. 



Taking the limit y — >■ cxj in this expression with the 
change of variables z — ky (and recalling that f{k = 
0) = 1), one deduces that in this limit: 



9{y) = - 



and consequently 



dk 



1 — cos(fcy) 



+ 0(l/y), (78) 
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dk (1 — cos(fcy)) 
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(79) 
Since the function (/(fc) — l)/fc^ is infinitely derivable, its 
Fourier transform vanishes in the limit y — >■ oo and one 
is left with 



m 



dk 



m - 1 



(80) 



(\ r°" dh r°° , o 

= -J 4 / dr sin{rk'/')rf{r) (81) 
TT Jo Vk Jo 
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dTT^/^fiT), (82) 
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using the definitions ([M)) . This universal quantity can 
thus be determined numerically with much higher pre- 
cision from F by performing a unique, non-oscillating 
integral. We find go ~ 1.149(18) which is in close agree- 
ment with the exact value, superior to the MC estimate 
go — 1.1137 [i^. All our estimates for the characteristic 
parameters of the different scaling functions obtained in 
this Section are summarized in Table HIl 



VII. CONCLUSION 

In this work, we have presented a general framework to 
investigate the KPZ equation using the NPRG method. 



We have proposed a detailed and revisited analysis of the 
symmetries of the KPZ equation, both the standard sym- 
metries - Galilean and shift - and their versions gauged 
in time, as well as the one-dimensional time-reversal sym- 
metry and a nonlinearly realized Z2 symmetry. We have 
derived general Ward identities associated with these 
symmetries (except for the Z2 one) and proposed a conve- 
nient geometric interpretation of the Galilean symmetry 
in terms of covariant time derivatives. 

We have then devised an approximation scheme based 
on an ansatz for r„ which, in contrast to the standard 
derivative expansion often used in the NPRG framework, 
preserves the momentum and frequency structure of the 
vertex functions. This ansatz, constrained by the sym- 
metries, allows one to compute correlation functions in 
a nonperturbative, yet systematically improvable, way. 
We explicitly implemented the minimal order in the re- 
sponse field of this approximation scheme, specifically 
developing two applications: the determination of criti- 
cal exponents in physical dimensions (with a simplified 
ansatz) and the computation of universal scaling func- 
tions in one dimension. 

We found that, without any input other than the bare 
action and its symmetries (in particular without assum- 
ing scaling), the renormalization group flow is generically 
attracted toward a strong-coupling fixed point - which 
roots the existence of generic scaling ~ or to the EW 
fixed point for d > 2 and A < Ac, yielding for the first 
time the full correct phase diagram of the KPZ equation 
within a RG approach. 

The estimates of the critical exponents obtained in di- 
mensions two and three compare reasonably with results 
from simulations, though the accuracy of these expo- 
nents decreases with the dimension. We are confident 
that these results, which were obtained using what re- 
mains, finally, a rather crude approximation, will get bet- 
ter as higher order approximations are studied - which 
is hopefully feasible in practice. Indeed, we know from 
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experience that within the NPRG framework, the con- 
vergence is generahy very fast, such that pushing it to 
even a shghtly higher order can eventuahy provide excel- 
lent results [30] . Similarly, we are unable, at this current 
order of approximation, to settle on the existence of an 
upper critical dimension, but we should be able to solve 
this issue at next orders. 

The scaling functions F, F and F obtained in one di- 
mension, on the other hand, compare very accurately 
with the exact ones. In particular, we recover the 
stretched exponential with tiny superimposed oscillations 
of the tail of f{k) on the correct scale fc^/^. We also ob- 
tained an accurate estimate of the universal amplitude 
ratio go ~ 1.149(18). Quantities beyond the leading scal- 
ing regime can also be computed, which we illustrated 
here with the determination of the correction to scaling 
exponent uj. 

A further step will consist of refining the computa- 
tion of critical exponents in all dimensions by improving 
the approximation, which is in progress. One could also 
compute the probability distributions of the h{t, x) field, 



which are again known exactly in one dimension but not 
in any other, and investigate the infiuence of boundary 
conditions. 
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Appendix A: Time-reversal Ward identities for the three- and four-point functions 



All the Ward identities for the three- and four-point vertex functions are derived by taking functional derivatives 
of the identity p7|) and evaluating them at uniform and static fields. For the remaining three-point function, one 
obtains in Fourier space: 

2r(f'°)(wi,pi;w2,P2) = -—plRerl^'^\uj2,P2;-^i -uj2,-pi -P2) - —plReT^i^'^\uJi,pi;-uji -uj2,-pi - P2) 



■P3R'er(f'^)(wi,pi;w2,_p2) 



I f v\^ 



PiP2P3rP^(^i:Pi;^2,P2). 



(83) 



For the four-point functions, we report the four additional independent identities obtained for uniform and static 
fields: 



2Imr(^2,2)(^^^^^^^^) ^ --^ppmrli'3)(^^^^2,^3) - -^p?Imr^3)(^^^^^^^^) 

- \(^)%lplplT^°'^\u^i.Lo^,u,2) - -^p^Reri2.2)(c,3,a;2,w4) - -^plReri2.2)(c,2, ,,4,0.3) 
-^pl'ReT^^^^\i^:i,uji,u:2) 



2Rer(,i'3)(^^^^2,W3) 
2Imrp)(cji,c^2,W3) 
2ReT'i'^\uj^,uj^,uj2) 



(84) 



and finally 



.(3,1) 



PilmT\. • •'(a;2,W3,-'^i - ^2 - ^3) + 3pcrm. = -- 



1 /l/\2 



2\D 



PiP2P3lnir^-^'^^(-i:Ji - W2 - 0Jz,i0i,L02) + 3 perm (85) 



where 'perm.' refers to the following: on the left-hand side, the permutation of the last (implicit) frequency- 
momentum argument (corresponding to the 1^ leg) together with the p^ factor in front, from wi over to 102, 1^3 and 
— cji — ijj2 — W3; on the right-hand side, the permutation of the first frequency- momentum argument (corresponding 
to the If leg) together with the complementary product of p^ in front, from — wi — uj2 — oj^ over to wi, uj2 and 0.3. 
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Appendix B: Validity of the approximation scheme an accurate description of the momentum dependence 



The approximation scheme presented in Section llVI is 
closely inspired by the BMW scheme, which provides 



of two-point functions at equilibrium [4l|. The BMW 
scheme thus stands as the most adequate approximation 
scheme for the KPZ problem, both because we are intrin- 
sically interested in the momentum and frequency struc- 
ture of the two-point functions and because the derivative 
nature of the interaction vertex seems to req uire one to 
go beyond a simple derivative expansion [39[. However, 
the implementation of the BMW scheme for the KPZ 
problem is hindered by the very demanding symmetries 
of the KPZ action. As a consequence, it cannot be di- 
rectly performed but has to be adapted regarding three 
main aspects. 

First, the standard BMW approximation at equilib- 
rium consists, at leading non trivial order, of perform- 
ing truncations on the three- and four- point functions, 
keeping only the leading order in their internal momenta 
since high momentum contributions are suppressed by 
the cutoff function (see ^41] for a full justification of 
the scheme). In the KPZ problem, as stressed in Sec- 
tion |lVl the Galilean symmetry relates n-point functions 
with different n and different momenta in a complicated 
way such that a naive expansion of the three- and four- 
point functions in their internal momenta and frequencies 
would spoil these relations. To overcome this difficulty, 
we propose here instead to devise an ansatz for P^ that 
manifestly satisfies all the symmetries and such that they 
are hence automatically preserved in subsequent calcula- 
tions (functional derivatives with respect to the fields). 
This ansatz imposes in turn the form of the three- and 
four- point functions. It appears that these functions re- 
tain some internal momentum and frequency dependen- 
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cies that would be neglected in a direct implementation 
of the BMW approximation but which are here necessary 
to fulfill the symmetry constraints. Hopefully, these ad- 
ditional subleading terms are not expected to deteriorate 
the quality of the approximation. 

Second, another specificity inherent to the KPZ prob- 
lem resides in the absence of regulator on the frequency 
sector, which is prevented by the Galilean-gauged sym- 
metry (see Section IIIID[) . As the whole BMW scheme 
relies on the presence of the cutoff term to perform expan- 
sions in internal momenta and frequencies, the absence 
of such a cutoff in the frequency sector would render an 
expansion in internal frequencies unjustified. Hopefully, 
the very same symmetry also cures this problem. Indeed, 
if vertices are expanded in momenta, then the Galilean- 
gauged symmetry fixes the frequency dependence associ- 
ated with ip-legs (see e.g. Eq. ([33])). Accordingly, no in- 
ternal frequency expansion is needed for the (/J-legs once 
an expansion in internal momenta is performed. 

Finally, the previous discussion applies to the (/s-legs. 
However, the (^-legs are treated in the proposed ansatz 
(|42p in a much cruder way since r„ is simply expanded 
at quadratic order in that field. The quality of such 
an expansion has been analyzed in equilibrium statisti- 
cal mechanics both within the derivative expansion (see, 
e.g. |2S,l3S]) and within the BMW approximation scheme 
[42| . and it seems to give reasonable results. Neverthe- 
less, it probably constitutes the main source of inaccu- 
racy of the present scheme. We leave for future work the 
improvement of the (p dependence of the approximation. 



Appendix C: n-point functions at the minimal order of the approximation scheme 



To compute the three- and four- point functions at vanishing field ensuing from the ansatz (|42|) involves to take 
functional derivatives of the functions f^ {—D^, — V^) {X — D, v, A) with respect to the field f and then to evaluate 
the result ai ip = 0. This can be done using the expression of the series expansion (j43p : 



m — l 






m— l.n— A;— 



(5(/3(ii,fi) 



(86) 



with 



= A {y^{5{t - ti)^('^)(f - fi)) • v,(at - Av,(p(x) • V,) 

+{dt - AV,^(x) • V,)(V,((5(i - ii)5('^)(f - fi)) • V,)) . (87) 



5ip{ti,xi 
When Fourier transforming. 



6ip{ti,xi) 



iXpi ■ {pi + p)(a;i + 2a;). 



One obtains for instance 



OO 771—1 

T'^y\u:i,pi,u:2,P2) = -i\ ^ a™™ x ^ ^f [^^ . ^^^ +j^^){uji + 20J2)] (wi + W2)'(""''"'H(pi +P2)')" + (2 ^ 3). 

m— l,n— k—0 

(89) 



and thus: 
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41-2) (c.i,pi, 6.2, P2) = -- [Pl • (Pl +P2) {f^iicOl + ^2)\ {Pl +P2?) - f^{u^l {pi +P2)')) 

-P,-P2 if^icolp^) f^{{uj,+u:2)^pi))] . 
All the other functions can be computed the same way. The three-point functions read 

r<;^'''HLOl,LU2;Pl,P2) = 

r^°^^\ivi,i02;pi.P2) - 

rl2.1)(wi,C^2;pi,p2) = Xp,.p2f^ {{lo, + u:2f, {Pl +P2f) + \pi.p2— [f^ {{^1 + UJ2)\pI) - f^ K,pt)] 

LO2 

+ *A^ U'^-(p*i +P2) [/: ((cji +CJ2)' , (Pl +P2f) - 1: [^i {pi+AY 



(90) 



(91) 

(92) 



-.2P2 -. 
Pl — -Pl 

W2 



!:({u:,+u:2f,pl)-i:{u:lpl) 



-.2P1 -. 

P2 ■P2 

Wl 



/:((a.i+c.2)',p|)~/:(. 



W2-P2) 



+ P2—-{Pl +P2) /k ((^1 + ^2)^ , (Pl +P2)^j - /^ (w|, (pi +P2)^j 



(93) 



(94) 



Note that the combinations appearing in Fk ' and Tk ' are related. If one denotes Fj;, ^ the expression (|90p where 

the D index labels the function (/^) appearing on the right-hand side, and similarly Fjj^i, ' the same expression 
where the function /^ is replaced by /^, one has 

iImr(f'^)(a;i,a;2;pi,P2) = -^ |pi^F(.^;j^)(w2, wi;P2,Pi) +p|r|^^;j^)(cji,cj2;Pi,P2)| • 

(31) 
This type of relation is used in the following to shorten the expression of F^ . The four-point functions are 

r'-^'°\uji,uj2,uj3;pi,P2,P3) = 

ri^'^^(wi,W2,W3;pi,P2,P3) = 

F(°'**)(a;i,a;2,a;3;pi,p2,P3) = 

r^^''^\uJi,UJ2,UJ3;Pl,P2,P3) = ~>?Pl-{Pl+P3)P2-PA 

-\^P2- {P2 + P3 )Pl -Pi 
-X^Pl. {pi + Pi )P2 -Pa 



I^iiiulpl) , /f (cji,Pl) ./k^ ((^1+^3)-, Pi 



(wi + a;2)w2 (wi + tJ2)wi 



a;iW2 



/f ((C.|,p|) , /f H,P|) /"i(-2+-3) ,pi 



^2 -2 



(wi -f 0)2)^1 (wi + 0^2)^2 



UJ1UJ2 



/f ((-I,p1) , /f (-ip|) ^" 1(^1 + "4) '^3 



\2 -2 



(wi + a;2)w2 (wi + tJ2)wi 



a;iW2 



-\^P2-{P2 +Pi)Pl-P3 



f^ii^i^pi) , f^H,pi) f^{i^2 + ^4)\pi) 



{UJI + UJ2)U}1 {UJI + CJ2)l^2 

— 7^(2,2)/ ^ ^ ^ ^ 

= i^K,D [^1^^2,^3;Pl,P2:P3) 



OJ1UJ2 



(95) 



where in the last line the D index refers to the involved function f^ . We denote F]^ ^' the two last lines of ((95)) with 
the same convention for the meaning of the index D. Finally, one finds, omitting the dependence in pi (which follows 
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the same order as Wi): 

rl3.1)(c,i,W2,C.3) = 



I 

UJ2 



^ -2r(2,2'M \ ^ -2T^i2.i)l \ ^ -2t^(2.2.)i \ 

-2J^P3^n:u'[^l,^2,^3)~ ^P2^K,u'[^3,i^l,i^2)- — Pi T), ,, ^ (cj2 , ^3, Wl) 

*^3f i^f ^ (^1,^2, W3) + i^2^^f^Jx^ (^3,^1,^2) + «Wif ^^f ^(CJ2, ^3, '^l) 

« P3-{Pl +P2)Pl-P2 ft {{^l +^2f , (Pl +P2)^j - Z;^ y^l, (Pl +P2)^j 

« Pl-(P2 +P3)P2-P3 /k ((W2 + Wa)^ , (P2 +P3)^j - Ik (^4) (P2 + Ps)^ j 

P2-(Pl +P3)Pl-P3 /k (('^l + ^3)^ , (Pl +P3)^j - Ik [^h (Pl +P3)j 

I 



(96) 



Appendix D: Procedure for the numerical 
integration 

This appendix is devoted to presenting the details of 
the numerical procedure implemented to achieve the in- 
tegration of flow equations, on the example of Eqs. ((52)) 
and ()53p . The momentum and frequency are discretized 
on a p X vw mesh of spacing Ap and A\/ro and sizes 
Pmax and V^max- Thc integrals over the internal mo- 
mentum q and frequency uj in (j5ip are performed using 
Simpson's rule. The vro^ mesh is chosen because the in- 
tegrand can have rather long tails in a) as there is no 
cutoff term (analogous to r{y)) suppressing thc high fre- 
quencies. The values of / at a) ± ro which do not fall on 
•\Aro mesh points are evaluated using cubic interpolations. 
The function / is extended outside the grid (for momenta 
p + q greater than Pmax and/or frequencies txi -I- a) greater 
than tiJmax) using power law extrapolations. 

Regarding the integration over q, as the integrand falls 
off exponentially due to the (derivative of the regulator) 
OsSk. term in ([?T|) . the bounds ±c» of the integral can 
be safely replaced by ±pmax- This is not the case for 
the integral over the frequency cu as the decay of the in- 
tegrand may be slow. We first compute the integral on 



[— tZJniax, 'f^max] and then evaluate the contribution of the 
integral on the boundaries [— oo, — ti7,nax] and [ti7niax,oo] 
by performing in these regions the change of variable 
X = Cj/Cj-ovuv. and using the extrapolated values of /. The 
precision on the momentum and frequency integral is of 
order 10~^ for the typical resolutions Ap — Avw = 1/4 
and mesh sizes w^ax = 15^ and Pmax from 20 to 45 in- 
creasing with the value of a (20 for a — 0.5 to 45 for 
a = 14). The derivative terms pdp and uid^ are com- 
puted using 5-point differences. 

We use an explicit Euler time stepping with a typical 
time step As — —2.10^^ to integrate the flow equations 
on the renormalization time s, which turns out to be sta- 
ble. In all cases, a fixed point is reached after s < —20. 
The fixed point functions are recorded at s = —25. We 
studied separately the influence of the resolution ( Ap and 
Avro^), and of the mesh sizes (pmax and tiimax) on the 
precision level, and checked the convergence. The dif- 
ferences (in all physical quantities) between resolutions 
and/or domain sizes are in all cases smaller than typi- 
cally 1% and these numerical errors are dominated by 
the residual variations observed when varying the cutoff 
parameter a. 
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